Observation of persistent centrosymmetricity in the hexagonal manganite family 
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The controversy regarding the ferroelectric behavior of hexagonal InMnOa is resolved by using a combination 
of x-ray diffraction (XRD), piezoresponse force microscopy (PFM), second harmonic generation (SHG), and 
density functional theory (DPT). While XRD data show a symmetry-lowering unit-cell tripling, which is also 
found in the multiferroic hexagonal manganites of P63 cm symmetry, PFM and SHG do not detect ferroelectricity 
at ambient or low temperature, in striking contrast to the behavior in the multiferroic counterparts. We propose 
instead a centro symmetric P3c phase as the ground state structure. Our DFT calculations reveal that the relative 
energy of the ferroelectric and nonferroelectric structures is determined by a competition between electrostatics 
and oxygen-/?- site covalency, with an absence of covalency favoring the ferroelectric phase. 



I. INTRODUCTION: STRUCTURE OF THE HEXAGONAL 
MANGANITES 

Hexagonal h-T^MnOs {R = Sc, Y, Dy-Lu) represents an es- 
tablished class of multiferroic s in which ferroelectricity and 
antiferromagnetism exist simultaneously. Although their fun- 
damental properties have been investigated for half a cen- 
tury, recent reports of intriguing characteristics such as in- 
terlocked antiphase (AP) and ferroelectric (FE) (AP+FE) 
domain walls^^ are fueling continued interest. At the 
root of these behaviors is their unusual improper geomet- 
ric ferroelectricity,^ ^ which is in turn related to their layered 
structure, in which xy planes of R^^ ions are interspaced by 
layers of corner-shared MnOs trigonal polyhedra (Fig.[T]). 

InMnOs crystallizes in the same hexagonal manganite 
structure as the h-7?Mn03 compounds, and might be expected 
to show analogous ferroelectric behavior. In fact, most previ- 
ous x-ray and neutron powder diffraction refinements assigned 
InMnOs to the polar PSscm structure adopted by the multi- 
ferroic hexagonal manganites.^ ^ In the PS^cm structure, the 
MnOs trigonal bipyramids tilt and trimerize with a trimeriza- 
tion phase of n • 60°, where n is an integer. The R ions on 
the 2a sites displace up or down along the z direction, depend- 
ing on the tilting direction, and those on the 4b sites in the 
opposite direction (Fig. [T]). This tilt symmetry then enables 
an additional displacement of the R sublattice relative to the 
Mn-0 layers causing a net ferroelectric polarization. Rusakov 
and Belik et al. pointed out that the nonpolar P3c structure - 
in which the MnOs polyhedra trimerize at intermediate an- 
gles and the inversion symmetry is retained (Fig. [T]) - and 
polar PS^cm structure have similar powder x-ray-diff'raction 
R values. 1^ They disregarded the P3c model in their subse- 
quent analysis, however, believing that all h-T^MnOs com- 
pounds should be polar. Indeed, ferroelectricity has been re- 
ported in InMnOs below 500 K based on the observation of 
polarization-electric field (P-E) hysteresis loops obtained by a 
ferroelectric test system. Such pyroelectric current measure- 
ments, especially if they are applied to amorphous samples or 
thin films, are notoriously sensitive to sample defects, how- 
ever, and the P-E loops shown in Ref. 10 could indicate leaky 
dielectric behavior rather than ferroelectricity.^^ In agreement 
with this, Belik et al. did not observe spontaneous polariza- 



tion when they repeated the experiment. Therefore there is 
no clear evidence to date that InMnOs has the PS^cm structure 
or shows ferroelectric polarization. 

In this study, we revisit the structure and polarization be- 
havior of InMnOs by using a combination of x-ray diff'rac- 
tion (XRD), piezoresponse force microscopy (PFM), optical 
second harmonic generation (SHG), and density functional 
theory (DFT) and show that InMnOs is indeed centrosym- 
metric, with P3c as the most likely space group (Fig. [l]). 
We explain the diff'erence between InMnOs and the multifer- 
roic h-7?Mn03 compounds using DFT analysis of the chemi- 
cal bonding, and propose another candidate material TlMnOs 
that should also show the nonferroelectric InMnOs structure. 
Finally, we discuss the implications of this newly-identified 
structure for the multiferroicity in the h-T^MnOs family in gen- 
eral. 



II. EXPERIMENTS 
A. XRD 

First we use powder XRD to directly and quantitatively 
compare the refinements for the candidate polar and nonpo- 
lar structures. For sample preparation, a stoichiometric mix- 
ture of In203 (99.9%) and Mn203 was placed in Au capsules 
and treated at 6 GPa in a belt- type high pressure apparatus 
at 1373 K for 30 min (heating rate 110 K/min). After heat 
treatment, the samples were quenched to room temperature, 
and the pressure was slowly released. The resultant samples 
were black dense pellets. Single-phase Mn203 was prepared 
from commercial Mn02 (99.99%) by heating in air at 923 K 
for 24 h. The synchrotron XRD data were obtained on pow- 
dered samples at the BL02B2 beamline of SPring-8.^^ They 
were collected in a 20 range from 2° to 75° with a step inter- 
val of 0.01 degrees and analyzed by the Rietveld method with 
RIETAN-2000.1^ 

An evaluation of the XRD data shown in Fig. [2] clearly re- 
veals that the unit cell of InMnOs holds six formula units. 
This indicates a deviation of the P63 /mmc high-temperature 
phase due to unit-cell tripling with tilt-shear motions of the 
MnOs bipyramids as in the ferroelectric T^MnOs compounds. 
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FIG. 1 : (a) Side and (b) top views of the two candidate InMnOs struc- 
tures. PG^cm is the ferroelectric phase of the h-/?Mn03 compounds; 
P3c is the structure proposed for InMnOa in this study. Arrows indi- 
cate the displacements from the high symmetry PG^/mmc phase. The 
primitive unit cell is shaded grey. In contrast with the PGscm phase, 
one In ion in the P3c unit cell remains at the high-symmetry 2b site, 
retaining the inversion symmetry. Note also the different tilt patterns 
of the MnOs polyhedra toward or around corner R ions, (c) Top 
views of Qf^ and /3~ domains in the PGscm phase using the notation 
from Ref. 1 and one of six domains in the P3c phase. The phases 0, 
defined by the (counter)clockwise angle of tilting direction of upper 
(lower) oxygen layers relative to a'^ domain, are also shown (Ref. 1 ). 
Note that the P3c phase has a tilt phase of 30°+f7-60°, and is obtained 
by averaging the tilt patterns and R-ion displacements of two PGscm 
trimerization domains with diff'erent origins and orientations such as 
the and /3~ domains. 



The centrosymmetric and the noncentrosymmetric subgroups 
with the highest possible symmetry that are compatible with 
a trimerization of the PS^/mmc high-temperature phase are 
PS^/mcm and PS^cm, respectively. Here refinements clearly 
favor the latter structure which may contribute to former 
claims of the PS^cm symmetry for InMnOs.^ However, in 
contrast to the case of the ferroelectric manganites, the struc- 
ture refinements of the XRD data reveals equally good fits 
for the centrosymmetric space group P3c and the noncen- 
trosymmetric space group PS^cm^^ consistent with a previ- 
ous observation.^ In Table |l| we report our refined atomic co- 
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FIG. 2: Synchrotron x-ray powder diffraction patterns of InMnOs 
at 293 K. Crosses represent data points and solid lines calculated 
intensities for P3c. Differences for three structure models are also 
shown. Bragg reflections are indicated by tick marks (these are the 
same for PG^cm, P3c, and PG^/mcm models). The lower tick marks 
indicate reflections from In203 impurity. Arrows show reflections 
corresponding to unit-cell tripling. 



TABLE I: Atomic fractional coordinates for InMnOa at 293 K 
in space group P3c (top) and P6^cm (bottom), obtained in this 
work using powder XRD. Our measured lattice constants a and c 
are 5.88462(10) and 11.48540(15) A for P3c and 5.88463(7) and 
1 1.48541(12) A for P6^cm, respectively. 
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FIG. 3: PFM images scaled relative to the response of periodically 
poled lithium niobate at room temperature, for (a) YMnOa show- 
ing the known characteristic ferroelectric domain pattern, and for (b) 
InMnOs, where no ferroelectric domains are observed, and only to- 
pographical artifacts (c) are visible. 



ordinates and lattice parameters within the P3c and PG^cm 
space groups. In general, structures refined with the correct 
space group are lower in energy than those refined with in- 
correct space groups. Our DFT calculations for InMnOs at 
the atomic positions and cell parameters obtained in the two 
competing best-fit experimental refinements (calculation de- 
tails given later) indicate that the nonpolar P3c structure is 
~200 meV per formula unit (f.u.) lower in energy than the 
PS^cm structure. We therefore suggest the nonferroelectric 
P3c phase as the ground state for InMnOs. 



B. PFM 

To confirm this suggestion, we next used PFM to probe di- 
rectly for the presence of FE domains; this technique avoids 
ambiguities caused by sample leakiness which might have oc- 
curred in previous macroscopic polarization measurements. 
In order to calibrate the response from InMnOs , PFM and si- 
multaneous scanning force microscopy (SFM) measurements 
with a commercial SFM (Solaris, NT-MDT) were carried out 
on YMnOs and InMnOs. All compounds for PFM and SHG 
measurements were grown by the flux method as z-oriented 
platelets. For PFM an ac-voltage of 14 Vpp at a frequency 
of 40 kHz was applied to a conductive Pt-Ir coated probe 
(NSC 35, Mikromasch). The out-of-plane component of the 
piezoelectric response was recorded by the in-phase output 
channel of an external lock-in amplifier (SR830, Stanford Re- 
search) with a typical sensitivity of 200 jjY and time constant 
of 10 ms. The PFM signal of each sample was normalized 
to a response of the z face of PPLN (^^33 =7 pm/V) in order to 
maintain comparability of the PFM response, which was mea- 
sured before and after each h-/?Mn03 measurement in order 
to exclude changes in the PFM sensitivity. 

Our results are summarized in the equally scaled Figs.|3ja) 
and [3jb). YMn03 reveals the familiar domain pattern of 
six intersecting AP-hFE domains with alternating polariza- 
tion ±pj^ Strikingly, the InMn03 shows an almost homo- 
geneous distribution of the PFM response with no sign of FE 
domains. The corresponding SFM data in Fig.[3jc) show that 
the slightly brighter speckles in the PFM image correspond to 




210 



180 



150 



210 



180 



150 



FIG. 4: SHG on (a) ErMnOg and (b) InMnOs at 5 K and 2.51 eV 
with light incident on the z-oriented surface of the single crystals 
along the [Oil] direction. The anisotropy measurement was obtained 
by rotating the linear polarization of the incident light at o) and of the 
detected intensity at loj simultaneously by 360°. Lines are fits of the 
SHG data according to Ref.[T8] Scales in (a) and (b) differ by a factor 
150. The two-fold SHG signal of ErMnOa indicates the presence of 
a spontaneous polarization whereas no indication for ferroelectricity 
in InMnOs was found. Instead the six-fold pattern characteristic of 
the antiferromagnetic order in InMnOs is observed. 



protrusions on the unpolished surface of the InMn03 sample. 
If one relates the contrast obtained for opposite domains in 
YMn03 to the spontaneous polarization of 5.6 yuC/cm^, any 
polarization in InMn03 has to be at least two orders of magni- 
tude smaller to avoid detection in our measurement. Our cal- 
culated polarization^^ for the DET-optimized PS^cm InMn03 
structure is 4.8 yuC/cm^ which would certainly be detectable. 
The absence of ferroelectricity (or of sub-resolution domains) 
in InMn03 is further supported by poling experiments with 
the SFM tip which did not induce any lasting change of the 
PFM response. 



C. SHG 

Since PFM measurements could only be done under ambi- 
ent conditions, we next used SHG to search for ferroelectric 
order at low temperature. As discussed in detail in Ref.W, the 
breaking of inversion symmetry by ferroelectric order leads 
to a characteristic SHG signal. The samples were mounted 
in a liquid-helium-operated cryostat and probed with 120 fs 
laser pulses in a standard transmission setup for SHG.^^ For 
comparison, ErMn03 was chosen for the SHG data because, 
unlike YMn03, it has the same magnetic SHG spectrum as 
InMn03 . Figure |4] shows the anisotropy of the SHG sig- 
nal taken under identical conditions at 5 K on ferroelectric 
ErMn03 and on InMn03. The laser light was incident un- 
der 45° to the hexagonal crystal axes so that SHG compo- 
nents coupling to a spontaneous polarization along z could be 
excited. ErMn03 shows the double lobe characteristic of the 
ferroelectric order. In InMn03 the double lobe is absent. In- 
stead a SHG signal with the sixfold anisotropy characteristic 
of the antiferromagnetic order and 150 times weaker intensity 
is found, indicating that the only order parameter is the anti- 
ferromagnetism of the Mn^^ ions. We therefore conclude that 
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FIG. 5: Calculated energy differences between the fully relaxed P3c 
phase and PGscm phases, A^str- Four exchange-correlation function- 
al, LDA, GGA, LDA+t/, and GGA+t/, and two spin configura- 
tions, FAFM and NCAFM (see text), were used. LDA/LDA+t/ re- 
sults for LuMnOa are not shown because vasp does not provide an 
LDA Lu PAW potential. Positive A£^str indicates that the ferroelec- 
tric PG^cm phase is stable over the P3c phase. The compounds are 
plotted in order of increasing R-sitc cation ionic radii, Rr (Ref. 20). 



InMnOi is not ferroelectric down to 5 K. 



III. THEORY 

To resolve the origin of the difference between InMnOs 
and the other hexagonal manganite h-T^MnOs compounds, 
we used DFT calculations to evaluate the energy difference, 
A£^str = Ep2c - Ep^^cm, between the candidate P6^cm and 
P?>c phases. Our spin-polarized first principles calculations 
were performed using the projector augmented- wave (PAW) 
method^ as implemented in vasp.^^ In this study, Sc 3^, ?>p, 
3d, and 4s, Y 4^, 4p, 4 J, and 5^, In 5^ and 5p, Lu 5p, 5d, 
and 6^', Tl 6^' and 6p, Mn 3d and 4^, and O 2s and 2p were 
described as valence electrons. The PAW data set with radial 
cutoffs of L3, 1.4, 1.6, 1.6, 1.7, 1.2, and 0.8 A, respectively, 
for Sc, Y, In, Lu, Tl, Mn and O was employed. The local 
density of states was also evaluated within the same spheres. 
Wave functions were expanded with plane waves up to an en- 
ergy cutoff of 500 eV. All calculations were performed with 
30-atom cells, which can describe unit cells of PS^cm and P3c 
phases, ^-points were sampled with a F-centered 4x4x2 grid. 
In addition to InMnOs, we also calculated A£^str for ScMnOs, 
LuMnOs, and YMnOs, as well as for as-yet-unsynthesized 
TlMnOs. To validate the results, we adopted four different 
exchange-correlation (XC) functionals: local density approx- 
imation (LDA), generalized gradient approximation (GGA), 
LDA-ht/, and GGA-hf/,^^ ^5 ^i^h the value for U^ff = U- J on 
the Mn-3 J orbitals set to 4 eV. In addition we tested two differ- 
ent spin configurations, so-called frustrated antiferromagnetic 
(FAFM)^^ and the noncollinear antiferromagnetic (NCAFM) 
adopted in Ref. 27 The lattice constants and internal positions 
were fully optimized in each case until the residual stresses 
and forces converged to less than 0.1 GPa and 0.01 eV/A re- 
spectively. 

Figure [5] shows our calculated A^str values for the various 
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FIG. 6: Calculated total and_averaged local densities of states for 
InMnOs and YMnOs in the P3c phase. Zeroes of the horizontal axes 
are set to the top of the valence band. Valence bands are composed 
mainly of 0-2p and Mn-3d orbitals. The difference of the electronic 
structures is mainly manifested at the formally unoccupied R5s and 
5p local density of states. 



functionals and magnetic configurations. A positive A£'str in- 
dicates that the ferroelectric structure is stable. When the R 
site is occupied with a IHb (Sc, Y, or Lu) ion, the PSscm phase 
is more stable than the P3c phase consistent with the experi- 
mentally observed ferroelectricity. However, in the case of a 
ma (In or Tl) ion, A£^str is close to zero. Note that in these 
calculations the lattice parameters and ionic positions for both 
structures were fully relaxed, within the constraint of the ap- 
propriate symmetry. [When we previously constrained our 
atomic positions and cell parameters to those obtained from 
the XRD analysis, the energy of the P3c phase is much lower 
(-200 meV/f.u.) than that of the PS^cm phase as we re- 



ported in Sec. |II A[ ] We note that A£^str is quite insensitive 
to the magnetic configuration, but does show a dependence 
on the choice of XC functional, and for InMnOs, the sign of 
A£^str depends on the functional, suggesting the possibility of 
competing low energy structures. [Previous calculations of 
the energy difference between P6^cm and P3c structures for 
InMnOs derived from XRD data (Sec.|II A[) and the polariza- 



tion for InMnOs with the P6^cm symmetry (Sec. II B), and 
subsequent calculations use the GGA-ht/ and FAFM configu- 
ration.] 

It is clear from Fig. [5] that Rr is not the key factor in de- 
termining the phase stability. Instead we focus on the differ- 
ent chemistry of the group ma ions, compared to the group 
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nb ions. It was previously suggested that the behavior of 
InMnOs is dominated by high-lying occupied semicore 4d 
(In) electrons in contrast the valence d states are formally 
unoccupied in the IHb ions. In fact it is well known that 
the presence or absence of semicore d electrons can affect 
the structural stability as illustrated by the different structures 
of MgO (rock salt) and ZnO (wurtzite, with semicore ds) in 
which the cations have very similar ionic radii (rMgo = 0-57 A 
and rzno = 0-60 A in four-coordination and rMgo = 0-72 A 
and rzno = 0-74 A in six-coordination). 

In Fig. |6ja) we show our calculated densities of 
states (DOS) for InMnOs and YMnOs (TlMnOs, and 
ScMnOs/LuMnOs behave analogously to InMnOs and 
YMnOs, respectively), both calculated within the P3c phase 
to allow a direct comparison. In both cases, the valence bands 
consist mainly of Mn-J (up-spin eig and e2g) and 0-2p states. 
The main differences occur in the DOSs on the R ions. The 
In "semicore" 4d states, however, form a narrow band that 
is around -13 eV below the top of the valence band when 
the d electrons are treated as valence (not shown). They do 
not directly contribute to covalent bonding with the oxygen 
anions, in contrast to the suggestion in Ref. 22, The rele- 
vant difference is the substantially lower energy of the for- 
mally unoccupied R 5s and 5p states in In compared with Y, 
caused by the well-known increase in nuclear charge without 
corresponding increase in screening across the 4d series. As 
a result, in InMnOs the 5^ (and to a lesser extent 5p) states, 
which would be completely empty in the ionic limit, develop 
significant occupation through ln-02p covalency, with occu- 
pied In 5 s states in fact forming the bottom of the valence 
band. (Si mila r behavior has been previously reported in other 
In oxides. ^^) In YMnOs, the Y 5s and 5p states are substan- 
tially higher in energy relative to the top of the valence band 
and so their hybridization with O 2p and subsequent occupa- 
tion is negligible. Instead there is a small hybridization with 
the formally empty Y 4 J states. 

The difference in covalency between InMnOs and YMnOs 
manifests particularly strikingly in the calculated valence 
charge densities at the R sites. In Fig.|7]we show the valence 
charge density differences between LuMnOs and YMnOs and 
between InMnOs and YMnOs in the P3c structure. The 
charge density at the Mn sites is similar in all cases. Com- 
pared with YMnOs and LuMnOs, however, InMnOs, has a 
decrease in charge density at the O sites adjacent to the In 
ions and an increase at the outer region of the In site indicat- 
ing charge transfer from oxygen to In and stronger In-0 than 
Y-0 or Lu-0 covalent bond formation. 

Next we investigate how the additional covalency of the In- 
O and Tl-0 bonds compared with those of Y-0 and related 
compounds manifest in the spring constants. Our calculated z- 
direction spring constants of the R ions at the high- symmetry 
2b sites in the P3c structure are 2.7, 3.2, and 3.1 eV/A^ for 
ScMnOs, LuMnOs, and YMnOs, and are 4.4 and 4.2 eV/A^ 
for InMnOs and TlMnOs . As expected, the strong In-0 and 
Tl-0 hybridization results in larger spring constants in the In 
and Tl compounds. The larger spring constants make the In 
and Tl ions reluctant to shift from their high- symmetry sites, 
favoring instead equal R-O^p bond distances. This in turn fa- 




FIG. 7: Charge density differences for valence band electrons be- 
tween (a) LuMnOs and YMnOs and (b) InMnOs and YMnOa in the 
P3c phase. For the comparison, the same atomic positions were cho- 
sen for the compounds that are compared. The structures are con- 
structed by averaging the DFT-optimized structures for (a) LuMnOs 
and YMnOa, and (b) InMnOa and YMnOs. Red and blue regions 
correspond to excess and deficiency of charge in LuMnOa/InMnOs 
compared with YMn03 . The yellow and light blue isosurfaces cor- 
respond to 0.03 and -0.03 A~^, respectively. Comparing with the 
charge density difference between LuMnOa and YMnOa, InMnOs 
shows noticeable differences from YMnOa especially at the In site. 



vors the P3c structure, in which | of the R ions retain their 
fully 6-coordinated high- symmetry positions, over the ferro- 
electric PS^cm phase, in which all R ions are displaced from 
the high- symmetry positions. 

We emphasize that the behavior here in which stronger co- 
valency favors the nonferroelectric phase is completely differ- 
ent from that in conventional ferroelectric s such as BaTiOs, 
in which stronger covalency favors the ferroelectric distor- 
tion through the second-order Jahn-Teller effect. In such con- 
ventional ferroelectric s, the Born effective charges Z* = ^ 
which participate actively in the re-hybridization are anoma- 
lously larger than the formal ionic charges, reflecting the 
charge transfer that takes place during the ionic displace- 
ments to the ferroelectric phase; such anomalous Born effec- 
tive charges are signatures of instability toward a ferroelectric 
phase transition.^ ^ In both InMnOs and YMnOs the mecha- 
nism for the primary symmetry-lowering tilt distortion is ge- 
ometric rather than due to a rehybridization, and the Born ef- 
fective charges on all atoms are nonanomalous.^^In InMnOs, 
the additional strong In-0 covalency in the paraelectric phase 
resists the distortion of the In ions away from their high- 
symmetry positions favoring the P3c space group, whereas 
the lower Y-0 covalency provides less resistance, allowing 
the additional Y-0 displacements required to reach the PS^cm 
symmetry. In Ref.|33l the hybridization between the Y-3d and 
0-2/7 orbitals was measured using polarization-dependent x- 
ray absorption spectroscopy (XAS) at the O ^-edge. Then the 
static charge occupancy in the Y 3d orbitals was equated with 
an anomalous dynamical Born effective charge, which led to 
the claim that this hybridization is responsible for the ferro- 
electricity in YMnOs. It is important to understand that the 
Born effective charge is the derivative of the polarization with 
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FIG. 8: Calculated (a) in-plane and (b) out-of-plane lattice constants 
in the P6^cm (o) and P3c (□) phases. The experimental lattice con- 
stants are also shown (x) (Ref. 6 34 35 ). The dashed lines indicate 
the trend of Htb manganites for a guide to eyes. 
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respect to ionic displacements, and is unrelated to the static 
orbital occupancy in a single structure: Partial hybridization 
of Y ?>d with 2/7, while clearly present both in the exper- 
iments and in earlier and subsequent first-principles calcula- 
tions, is not indicative of an anomalous Born eff'ective charge 
and therefore does not indicate tendency toward ferroelectric- 
ity. 

Experimentally, it is known that InMnOs has an anoma- 
lously large c lattice constant compared to the multiferroic 
hexagonal manganites.^ In Fig. [s] we plot our calculated lat- 
tice constants (with experimental values where available for 
comparison) of the manganites series as a function of ionic 
radii Rr. We point out first that this is not a consequence of 
the diff'erent space group that we have established here; our 
density functional calculations yield similar lattice constants 
for InMnOs in the P?>c and PS^cm phases. The calculated 
lattice constants are systematically overestimated compared 
with experiments as is typical of the GGA. We see that the 
in-plane lattice constants, a, increase monotonically with Rr, 
with InMnOs showing only a small calculated anomaly. In 
contrast, the c lattice constant of InMnOs deviates strongly 
from the trend shown by the IKb manganites, both in our cal- 
culations and in experiment. Since this deviation is also iden- 
tified in TlMnOs, the anomalously large c likely originates 
from covalency in IHa manganites as shown in Fig.|7] 

IV. DISCUSSION 

Although the centrosymmetric P3c phase that we propose 
in this work for InMnOs may be seemingly less attractive 
compared with the ferroelectric PS^cm structure, our results 
have implications for the multiferroic hexagonal manganites 
as a whole. Since the tilt pattern of the YMnOs structure 
subsequently allows for the development of ferroelectricity, 
whereas that of the InMnOs structure does not, the subtle 
chemical bonding diff'erences identified here that favor one tilt 
pattern over another in turn determine whether the resulting 
structure can be multiferroic. Specifically, we have discussed 
here that an absence of R-0 hybridization is required to favor 
the YMnOs tilt pattern over the InMnOs tilt pattern; an ab- 



FIG. 9: (a),(b) Calculated InMnOa P3c and P6^cm structure energies 
relative to the equilibrium P3c phase as a function of in-plane lat- 
tice constant a by the^GGA and GGA+t/. (c),(d) Energy differences 
A^str between the P3c and PG^cm phases. Positive values indicate 
that the ferroelectric P6^cm phase is stable over the P3c phase. The 
fitted lines were obtained with fourth-order polynomials. The verti- 
cal dashed lines indicate the calculated equilibrium lattice constants 
of InMnOs P3c phase. 



sence of R-0 hybridization is therefore a requirement for fer- 
roelectricity in the hexagonal manganites. Earlier theoretical 
papers correctly noted the electrostatic origin of the "geomet- 
ric ferroelectricity" mechanism in the hexagonal manganites 
we now understand that the relative stability of ferroelectric 
and nonferroelectric structures is determined by a competi- 
tion between electrostatics (favoring the ferroelectric phase) 
and covalency (favoring the nonferroelectric phase). 

In addition, since the polar PS^cm and nonpolar P3c phases 
are close in energy in InMnOs as shown in Fig. [5] it might be 
expected that their relative stability could be changed using 
external perturbations such as epitaxial strain. Figure |9] shows 
the calculated energies of the PSscm and P3c phases and their 
energy diff'erences as a function of in-plane lattice constant a, 
with the out-of-plane lattice constant c and internal positions 
fully relaxed for each a value. Because the critical strain of 
the phase boundary depends on the U value, we performed 
the calculations using both the GGA and GGA+f/ methods. 
We see that the A^str increases with the increasing in-plane 
lattice constant, indicating a larger in-plane lattice constant 
could develop the polar P63 cm. Interestingly, the ferroelectric 
polarization develops in the out-of-plane direction, in striking 
contrast to the behavior in perovskites.^^ Therefore we antic- 
ipate that InMnOs could be tuned into the polar PS^cm struc- 
ture using tensile strain. 

Finally we mention that a recent transmission electron mi- 
croscopy study of the domain walls in ferroelectric hexag- 
onal TmMnOs and LuMnOs^^ revealed that a domain wall 
structure at the edges of the sample is similar to the nonpolar 
P3c InMnOs structure: The R ion at the wall is at the cen- 
trosymmetric position, with one neighbor displaced in the up- 
direction and one in the down-direction. Detailed calculations 
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of the domain- wall structure in /^MnOs are ongoing. 

V. CONCLUSION 

In summary, we have proposed a different nonferroelectric 
ground state structure in the hexagonal manganite InMnOs, 
and we predict its occurrence in as-yet-unsynthesized hexag- 
onal TlMnOs. The proposed phase has P3c symmetry, and is 
closely related to the usual P63 cm ferroelectric ground state 
but with a different pattern of polyhedral tilts that retains the 
center of inversion. The energy balance between the two re- 
lated phases is determined by a competition between electro- 
statics and R-0 co valency, with lower R-0 covalency favoring 
the ferroelectric structure. Thus, the absence of ferroelectric- 
ity in InMnOs reveals to us the reason for the presence of 



ferroelectricity (and therefore multiferroicity) in the other h- 
T^MnOs compounds. 
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